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Abstract 

In this paper we consider the bilocal field approach for QCD. We obtain 
a bilocal effective meson action with a potential kernel given in relativistic 
covariant form. The corresponding Schwinger-Dyson and Bethe-Salpeter 
equations are investigated in detail. By introducing weak interactions into 
the theory we study heavy meson properties as decay constants and semilep- 
tonic decay amplitudes. Thereby, the transition from the bilocal field de- 
scription to the heavy quark effective theory is discussed. Considering as 
example the semileptonic decay of a pseudoscalar B-meson into a pseu- 
doscalar .D-meson we obtain an integral expression for the corresponding 
Isgur-Wise function in terms of meson wave functions. 
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1. Introduction 

The investigation of heavy meson decays has become one of the important problems 
in heavy flavour physics. Especially, -B-decays [1] play an important role in determining 
the Kobayashi-Maskawa matrix elements, including the CP violating phase. Further- 
more, rare decays of heavy mesons may indicate deviations from the standard model. 
For the description of physical (hadronic) decay processes one needs to know the wave 
functions and form factors. There exist several approaches to attack these problems. 
Examples are phenomenological quark models, QCD sum rules, heavy quark effective 
theory, potential models. 

Here we shall consider one more possibility - potential models in the bilocal field 
approach for QCD [2,3]. Thereby, we start from an approximate QCD action with 
massive quark fields and hadronize it. Substitution of the interaction kernel containing 
the free gluon propagator by a instantaneous Lorentz-vector potential kernel yields a 
class of bilocal potential models for meson fields [4, 5]. Potential models do not possess 
full relativistic invariance because the Fock space is restricted to (qq) pairs and an 
instantaneous interaction is assumed. Nevertheless, in our model the kernel is written 
in relativistic covariant form by introducing a special vector being proportional to the 
bound-state total momentum. This allows us to investigate equations for a bound state 
moving together with the potential kernel. This fact is important especially for the 
calculation of formfactors where one needs to know meson wave functions for moving 
particles. One further advantage of our relativistic model consists in the possibility 
of describing in a uniform manner both light and heavy mesons in dependence on 
the choice of the potential. So, in the case of a short-range potential one can use 
the separable approximation (for small orbital momenta) to obtain a new regularized 
version of the Nambu-Jona-Lasinio model [6]. We should also mention that we are able 
to extend some of the results obtained by A. Le Yaouanc et al. [7], who investigated 
a local quark model for massless quark fields within the Hamilton formalism with an 
instantaneous fourth-component Lorentz-vector colour confining potential (e.g. the 
harmonic oscillator potential). 

In the first part of this paper the Schwinger-Dyson and Bethe-Salpeter equations 
for our model will be derived, and the meson wave functions fulfilling Schrodinger-like 
equations will be introduced. Thereby, the equations for the quark mass spectrum will 
be considered in application to oscillator, Coulomb and linear potentials. Concerning 
the equations on bound state masses special attention will be paid to the derivation of 
equations for wave functions of scalar, pseudoscalar and vector particles. 

Then, in the second part of the paper we will apply the bilocal field method to heavy 
meson physics. Therefore we will introduce semileptonic weak interactions into the 
underlying potential theory by shifting the bilocal field correspondingly. Furthermore, 
now we will have to consider in general moving bound states. We will derive formulas 
for pseudoscalar meson decay constants as well as for semileptonic decays of heavy 
quarkonia. Moreover, we will establish a relation between bilocal field approach and 
heavy quark effective theory. This will allow us to consider the semileptonic decays 
in the limit to the heavy quark effective theory [8 — 10]. Thereby, as example the 
Isgur-Wise function [8] corresponding to decays B — > D(lvi) will be defined within the 
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bilocal field method. 

In this paper general formulas containing the so far undetermined meson wave 
functions will be obtained. Nevertheless, as will be shown, these wave functions fulfil 
systems of integral equations which shall be solved in the near future for concrete 
potentials by different methods. 

The paper is organized as follows. In sect. 2 the formulation of the model is given. 
The corresponding equations for the quark spectrum within a meson are derived in 
sect. 3. Sect. 4 contains the investigation of equations for bound state vertex and 
wave functions. In sect. 5 formulas for the pseudoscalar meson decay constants are 
derived. In sect. 6 we discuss the relation between bilocal field approach and heavy 
quark effective theory. Semileptonic decays of heavy quarkonia in the limit of heavy 
quark effective theory are investigated in sect. 7. The summary and conclusions are 
given in sect. 8. 

2. Formulation of the model 

2.1. Hadronization of QCD 

Let us start with the approximate QCD action for quarks q in the form 

2 

W = J d 4 xq(x)[G \x)]q(x) - 9 - j j d 4 xd 4 yj ;(x)Df u (x - y)j b v (y) . (1) 

Here G is the Green's function for free quarks, 

Gq 1 = ip -rh , 

where rho is the bare quark mass matrix, rh = diag(m°, m!], . . . ,m% ), Nf being the 
flavour number. The quark current j£(x) is defined by the relation 

where A a are the Gell-Mann matrices in colour space SU(3) , // is the Lorentz index 
and 7^ denotes the Dirac matrix. The quark-gluon interaction with coupling constant 
g is mediated via the free gluon propagator D^(x — y) = 5 ab g^ u D(x — y). 

For hadronization of action (1) let us first rewrite the bilocal interaction term 

- J j ' d*xdS\(q{x)^q{x)) 9 ^D{x - y)(q(y)^q(y)) 

in the form 

d 4 xd i yq B {y)qA(x)K ABtEF (x - y)qF(x)q E (y) (2) 
with the kernel 

KabM* -V) = ^Mts E ^f^f 5uS kj 9 ^D(x - y) . 
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Here A, B, E, F are short-hand notations for the indices A = {r, a, i}, B = {s, /3,j}, 
E = {£, 7, k} and F = {u, 5, 1} . The first index in the bracket refers to the Lorentz 
group , the second one to the colour group SU(3) C and the third one to the flavour 
group SU(N f )f. 

Now, we make the colour Fierz rearrangement [11] 

8 ^ 2 3 

H KsKp = 0^/^7*5 + oE e pa~f£ P ^ (3) 
o=l 6 6 p=l 

with e Q/ 3 7 being the antisymmetric Levi- Civita tensor. This identity allows to de- 
compose Kab,ef(% — y) completely into "attractive" colour singlet (qq) and colour 
antitriplet (qq) channels. In this treatment "repulsive" colour octet (qq) and sextet 
(qq) channels are absent in a natural way. In this paper we want to discuss (qq) meson 
bound states. For this reason we will consider only the colour singlet part of (3) and 
therefore the kernel 

■ i g 2 

k ab 9 ef(x ~y) = lru(l^tsj^S a ^ 5 5ii5 kj —D(x - y) , (4) 

where N c = 3 denotes the colour number. 

Inserting the kernel (4) into the bilocal interaction term (2) one obtains (omitting 
the group indices) 

w 1 1 ^^yiiy^h^ ^ ~ y^ q ( x>) y(y^ ■ 

Then, action (1) can be represented in the form 

W M = J jd*xd 4 y{(q(y)q(x))(-G^)5(x-y) 

+ 2~k [( 9 (^(aO)£(s - y)(q(x)q(y)j\ } 

with the Lorentz- vector kernel for QCD 

K,(x-y)=g 2 1 »D(x-y)® lll . (5) 
In symbolic notation it reads 

W M = (qq,-G^) + ^(qq,JCqq) . 

2.2. Schwinger-Dyson and Bethe-Salpeter equations 

Let us now consider the functional integral 

Z = J VqVqexp{iW M [q,q}} . 
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After integrating over the quark fields with the help of the Legendre transform one 
gets 



Z = J VMexp{iW eff [M}} , 



with the effective action 



W eff [M] = N^-^M^-'M) - tTM-G, 1 + M)} (6) 

for meson fields Ai. Here Tr means integration over the continuous variables and 
taking the traces over the discret ones (spinor and flavour indices). The condition of 
minimum for this effective action reads 

Let us denote the solution of this equation by (E — rh°). Then one obtains from (7) 
the Schwinger-Dyson equation 

E = m° + i/CG s , (8) 

where 

Gs 1 = ip - E . 

Expanding action (6) around the minimum with M. = (E — m°) + M! one gets 

w eff [M] = w eff m 

r 1 7 00 1 i 

+ N e \ --(M', IC-'M') - -Tr(GxM') 2 - i Y, -Tr(-G s A^') n ■ (9) 
I z Z n=3 n ) 

The vanishing of the second variation of this effective action with respect to M.' , 



S 2 W P 



• r = o, 

M'=0 



5M'5M' 

leads to the Bethe-Salpeter equation 

T = -iKiGvTG*) (10) 

for the vertex function F(x,y) in the ladder approximation. 

The Schwinger-Dyson equation (8) describes the quark spectrum in the meson, 
whereas the Bethe-Salpeter equation (10) yields the bound state spectrum. Solving 
together both equations one may obtain the wave functions of the bound states and 
calculate with their help not only static properties of mesons as the mass spectrum 
and the decay constants but also decay probabilities. 
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2.3. Relativistic covariant description of potential models 

From the translation invariance of the two-particle bound state it follows that one 
can separate the cm. motion of the (qq) system from the relative motion. But in 
a relativistic theory it is impossible to separate the cm. coordinates. Nevertheless, 
one may separate the total momentum being the momentum of the cm. motion. 
After this one can more or less arbitrarily define some coordinate representing the 
absolute position in space-time and use the relative coordinate of the bound state only. 
For example, let X^ be the position x or y of one of the particles (quarks) or a linear 
combination X = ax+ (l — a)y of them. For a — 1 the relative coordinate is z = x — y. 
In this case one can look for the solution of the Bethe-Salpeter equation for the bound 
state wave function ip in the form 

Let us now substitute the QCD kernel K, of the integral equations (8) and (10) by 
a potential one. This can be done in a relativistic covariant way by choosing instead 
of (5) 

lC(x, y) = K?(x - y\^-) = iV(z L )5(z ■ v ) <g> fi (11) 

with 

Zp=z»-zl ' z l = ^vk z ■ ■ 

Here f]^ is a vector (fj = 77^7^; rf = 1), being proportional to the momentum eigenvector 
and describing the motion of the bound state as a whole 

In (11) the transversality of the exchange interaction in the (qq) system is ensured 
by V(z ± ) - some phenomenological potential for the description of quarkonia. Fur- 
thermore, the S— function S(z ■ rj) guarantees the instantaneousness of the exchange 
interaction. 

We should add that one can arrive at equation (11) by discussing moving bound 
states rigorously within the quantization theory for gauge theories [12]. For the bound 
state at rest one has V = 0, so that rj^ = (1, 0, 0, 0) and the kernel takes the form 

/C = 7oVV)<H^o)®7o • ( 12 ) 

Kernel (12) is known from the calculation of the (e + e~) positronium spectrum in elec- 
trodynamics with V as Coulomb potential, where an exchange of transversal photons 
takes place. 

Phenomenologically we may choose the interaction potential in the form 

4 a 

V(r) = --— + V 5( y r) + ar + br 2 . (13) 
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Here the first term is a Coulomb-type potential for one-gluon exchange, whereas the 
last two terms in (13) guarantee quark confinement. The second term leads to the 
Nambu-Jona-Lasinio model which has been considered within this approach in [6]. 

The potential (13) may be attributed to the sum of the Coulomb and oscillator 
potentials. In the applications we will consider only two potentials V(r), the first one 
being the sum of Coulomb, constant and linear potentials and the second one - the 
sum of Coulomb, constant and oscillator potentials. In Table 1. all potentials under 
consideration are displayed in x space as well as in momentum space. (There the 
definition A p = d 2 /dp 2 has been used.) 

3. Equations for the quark mass spectrum within a meson 

3.1. The case of an unspecified potential 

The Schwinger-Dyson equation (8) takes for the potential kernel (11) in momentum 
space the form 



Assuming flavour diagonality of S(p), i.e. considering £ = diag(Ei, E2, ^N f ) , 
equation (14) falls into identical equations for all S n with bare mass term mJJ , n = 
\,2,...,Nf. Omitting the flavour index n, these Schwinger-Dyson equations for the 
fields E of a given flavour read in the rest frame, r/^ = (1, 0, 0, 0): 




(14) 



where 




S(P) 



m + i 




'0 



(15) 



Now we make the ansatz 




1 = 



1,2,3. 



(16) 



Then, the propagator Gs(g) for the free quark can be represented as follows 



1 



A+(q) A_(q) 



7o 



^-£(q)+*e 



g - £(q) + ie go + E(q) - ie 



(17) 
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Here the notations 



A±(q) = ^ 1 (q)A° ± 5(q) = ^(l±^ 2 (q)7o) = i(l±7o5 2 (q)) 



A ± (q) = S( q )A ± S-\ q ) = 1(1 ± 5 2 (q) 7o ) = ^(1 ± loS-\q)) 



(18) 



and 



^(q) = |q|[A 2 (q) + (l +J B(q)) 2 ] 1/2 
have been introduced, where 

S ±2 (q) = sin0(q) ± qcos0(q) = exp{±2qz/(q)} , 



(19) 



sin0(q) 



^(q)|q| 



cos0(q) 



(1 + S(q))|q| 



£(q) ' ™^ v ^ E(q) 
S^q) = cosz/(q) ± qsinz/(q) , z/(q) = ±(-0(q) + f ) , 



q = iiii , ii = 



-l 



(20) 



and 



A 



±-^(l±7o) 



(21) 



With these definitions one can write the solution S(p) of the Schwinger-Dyson 
equation (15) in the form 

E(p) = £(p)sin0(p) + p(£(p)cos0(p) - |p|) = (p 7 ) + E(p)S- 2 (p) . (22) 

By inserting (17) and (22) into (15) we can rewrite the Schwinger-Dyson equation as 

£(p)sin0(p) - p(£(p)cos0(p) - |p|) 

= rn° - - J J^yV(P ~ q)7o(l ~ 7o(sin0(q) + qcos0(q))) . 

Taking the trace on both sides of this equation, one obtains a system of two equations, 



£'(p)sin0(p) = m° 



rfq 



V(p - q)sin0(q) 



2 J (2tt) 3 
1 f do 

E(p)cos0(p) = |p| + - J 77^^(P - q)x(P, q)cos0(q) , 



(23) 



which defines the mass spectrum of two quarks forming a bound state. Here the 
notation 



x(P,q) = (q- P) = cos#(p,q) 



(24) 



has been introduced. For solving the system of equations (23) it is necessary to fix the 
interaction potential. 

3.2. Application to concrete potentials 

Let us now investigate the system of equations (23) obtained from the Schwinger- 
Dyson equation (15) for different potentials. First of all we rewrite the potential (13) 
in momentum space: 



4 Att Rir 

V(P - q) = -i^s)j—y - «k=^4 + V ° ~ K2vr) 3 A q 5( P - q) . (25) 



The last term in this potential - the oscillator term - should be considered separately. 
Inserting it into the system (23) one gets 



b f 

£(p)sin0( P ) = m - - J dq(A q 5(p - q))sin0(q) 

b r 

^ £(p)cos0(p) = |p| - g / rfq(A q 5(p - q))x(p, q)cos0(q) 



(26) 



Now we make use of the relations 

2 

A q (sin0(q))| p=q = -cos0 • 0' - sin0 • (0') 2 + cos0 • (0") , 

2 2 

A q (cos6'( P , q) • sin0(q))| p=q = cos^ • ( -coscp sin0 • (0 ) 

-cos0 • (0') 2 - sin0(0")) , 

where the notation p = |p| has been used. Then, the system (26) corresponding to the 
Schwinger-Dyson equation for the oscillator potential leads to the differential equation 

b . b 

— -(p 2 <p r )' = p 3 sin0 — mop 2 cos0 + -sin(20) 

on the function 0(p).Q Knowing it's solution one may in principle calculate the con- 
stituent quark mass m in dependence on two parameters - the current mass m° and 
the potential parameter b. Indeed, let us consider relations (16), (19), (20) for the case 
A(p) = m/\p\, -B(p) = , so that one has 



£(p) = m, E(p) = yp 2 + m? , 

(27) 

cos0(p) = — j- , sin0(p) 



E(p) ' — E(p) 

Here the last relation is the equation on the constituent mass m. 



^^For mo = this equation has already been given in [7]. 
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Now we turn to the linear and Coulomb potentials in (25). In this case one has 
to investigate the renormalized Schwinger-Dyson equation because the second equa- 
tion of system (23) contains an ultraviolet divergence. We introduce into the latter a 
renormalization constant Z, 

1 f dq 

£(p)cos</>(p) = Z\p\ + -J 77^V(P - q) X(P, q) cos0(q) . 



and choose it as 



1 - 



(2k) 
dq 



/?2^ (p - q)x(p ' q) - 



2|p| J (2tt) 3 

Then the system (23) can be represented in the following manner 

1 f dq 

£(p)sin0(p) =mQ + -J __y(p _ q )sin0(q) 

* 1 f dq 

£(p)cos0(p) = |p| + - j J^ V (P ~ q)x(P, q)(cos0(g) - 1) . 

From here the energy E(p) may also be expressed via the function 0(p): 
E(p) 



(28) 



+ 



m o sin0(p) + |p|cos0(p) 
dq 



V(p - q){sin0(p)sin0(q) - X (p, q)cos0(p)(cos0(q) - 1)} . (29) 



The infrared singularity appearing in (28) and (29) for the Coulomb potential may 
be removed if one changes the physical observable - the excitation energy AE(p) = 
E(p)-E(0). 

Furthermore, from (28) one obtains the renormalized integral equation 



|p|sin0(p) = m + 



1 



dq 



V(p — q){cos0(p)sin0(q) 



-X(p,q)sin0(p)(cos0(q) - 1)} 



(30) 



Now, we integrate in (30) over the angular variables. Then, for the Coulomb potential 
one gets 



|p|sin0(p) 



m o cos0(p) + 



1 4a„ 



2?r 3 

is! + 1 q 2 + p 2 



d\q\ 



In 



|q| - 


IpI 


|q| + 


IpI 



cos0(p)sin0(q) 



In 



|q| - 


IpI 


|q| + 


IpI 



sin0(p)(cos0(q) — 1) > . 



In the case of the linear potential equation (30) becomes 

|p|sin0(p) = m o cos0(p) - J rf|q|| ^ q cos<ft(p)sin0(q) 

I0(p)cos0(q)| 



|q| q 2 + p 2 i. 



IpI (q 2 -p 2 ) 2 2 P 2 



|q| - 


IpI 


|q| + 


IpI 



10 



The singularity in this equation can be regularized by introducing a small cut-off A 
around the singular point p = q and taking the limit A — > afterwards. 

4. Equations on bound state masses 

4.1. Equations for bound state vertex functions 

The Bethe-Salpeter equation (10) for the vertex functions T reads in the case of 
the potential kernel (11) in momentum space 

r( P |P) = -iJ^Vip-^G^q + ^y^G^q-^f,, (31) 

where G n = G^ n is given by (17). Here the index n = 1, 2 is used to distinguish between 
the two quarks forming the bound state. Therefore the quantities E, A, B, A, S, E, 0, v 
defined by (15)-(20) will now carry this index. The quantity V denotes as before the 
total momentum of the bound state. 

For what follows we will investigate equation (31) in the rest frame, 77^ = (1, 0, 0, 0): 

t(p) = -if t|^(p - qhoG^g + f )r( q )G 2 ( g - ^) 7o . (32) 

Here M means the bound state mass. Then, according to (17) the Green's functions 
G n may be expressed as 

G(a ±M = ( A + } (q) | A - } (q) V. 

n{q 2 1 \q ±M/2-E n (ci)+ie q ± M/2 + E n (q) - it) 10 

(33) 

Af)(q) A W (q) \ 



= 7o I 



go ± M/2 - £ n (q) + te q ± M/2 + E n (q) -iej ' 



whereby E n (q),n =1,2 are the solutions of the system of equations (23). Integrating 
the Bethe-Salpeter equation (32) over go we obtain the Salpeter equation 

with E(q) = £?i(q) + E 2 (q) and 

n ±T ( q )=Ai 1) ( q ) 7 or 7 oA? ) (q) • (35) 
Let us now decompose T(p) with respect to it's Dirac structure, 

r = r 1 + 7o r 2 , 

where 

r, = 7 5 -rf + 7 p -rf + 7 rrr + 7^ , / = i, 2 , < = i,2,3. 



n 



Here we have defined 7 5 = 1, 7 P = 75, 7^ 



7*> 7-, 



,4 



7i75 for scalar, pseudoscalar, 



vector and axial-vector bound states T s , T p , T v , T A , respectively. It is favourable to 
write T in the form 



(36) 



so that the expression (7 r7 ) in (34) may be rewritten as 



with 



( 7o r 7 o) = E «/(r( + 7oi1)7 / 

J=S,P,V,A 



{«/, / = 5,p,y,A} = {i,-i,-i,i} 



(37) 



Then, inserting (37) into (35) one can derive from (34) the Bethe-Salpeter equations 
for the vertex functions r[. Here we shall give only the result of the calculation for 
the simple case, in which the relations (27) hold for both quarks. Then the Green's 
functions (17) are given by G n (q) = — m n + ie)~ l , n — 1, 2 with m n as constituent 
quark mass. One obtains 



r£(p) 



+ (1 - a//5 7 )q 2 



E\ Ei 



r((q) + M 



mi m 2 
Ei Ei J 



+ 



nP-c)^{M(|i-^)r;( q ) 

(mi - «^)(| -§) + (! + «'A)1 2 (^ + ^)] r =(l) 



with £ n = £„(q) 



' + , n = 1, 2, E = Ei + E 2 and 

{ft, J = s,P,v,A} = {-i,i,i —} . 



4.2. Equations for bound state wave functions 

Sometimes it is favourable to work not with the vertex function T but with the 
Bethe-Salpeter wave function For our potential theory (9), (11) both quantities are 
connected with each other in the rest frame by the relation 

r(p) = / (2^ (p " q ) 7o ^ q ) 70 ' ( 38 ) 
such that according to (34) ^ is defined as 

¥(q) = n+ - (q) + n - +(q) (39) 
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with n± T (q) given by (35). Furthermore, it is more suitable to introduce instaed of 
(39) a new wave function 

I (q) = Si(q)tt(q)S 2 (q) ■ (40) 

Then one may express the quantities H± T (q) appearing in the definition (39) of the 
wave function ^(q) as follows: 

n ±T (q) = S^(d) n ±T (q)S 2 \ q ) (41) 



with 



o ooo 

n ±T (q) =A±r (q) A T , (42) 



o 

where A± is defined by (21) and 

f (q) = 5r 1 (q)r(q)5 2 - 1 (q) . (43) 

With the help of these relations one obtains from (38) a Salpeter equation for the new 
quantities. It reads 

f (p) = - / (0^(p - <i)^(p> q) * (q)^(q, p) , (44) 

where 

S;(p,q) = 5- 1 (p) 7o 5- 1 (q), n=l,2. 

Notice, that according to (40) and (41) the wave function (q) has a representation 
in analogy to relation (39) for ^(q): 

* (q) - j + ; { t + L+(q) 



E(q) - M E(q) + M 
From here one gets the two relations 



(q) l ^ E(^M ^ (q) ' (45) 
Let us rewrite them in a form similar to Schrodinger equations 

[E(q) T M] A±4 (q) A T =n ±T (q) . (46) 
o 

Next, we decompose the wave function ^ (q) of two-particle bound states, 

oo o 
^=\Pi +7o- *2 , 

13 







u 

and expand I — 1, 2 in analogy to (36) over the full system {7 J , J = 1,2,3} = 
{l5,e a (p),p}, so that 
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tf=£[*i+7o-*£]7 J - (47) 
j=i 

After inserting decomposition (47) into the Schrodinger-like equations (46) one may 

o o 

derive the following systems of equations on and \I/ 2 f° r ^ = 1> 2, 3: 
M k (p) • itr( 7 V) = ^(P) *i (P) • |tr( 7 V) 

M ^ (p) • ltr( 7 V) = #(p) k (P) ' ^tr( 7 V) 



(48) 



where 



^ L (p,q) = ^r( 7 ^^(p,q)7 L ^(q,P)) 



The calculation of the traces yields 

T£ L (p,q) = c^c^tr( 7 L 7 X ) - s^^tr( 7 L 7 x pq) 

with 

{ PK \K= 1,2, 3} = {-1,-1,1} 

and 

c% pK = c 2 (q)ci(q) =FPA-S2(q)si(q) , 
= s 2(q)ci(q) TPi^s 2 (q)ci(q) , 
s„(q) = sin0 n (q) , c„(q) = cos0„(q) , n = 1,2 . 

o 

The systems (48) of integral equations for the bound state wave functions v&f and 

o 

#2 nee d to be specified for every particle type. To do this let us introduce the notations 

o o 

*]=Li, , *i=Vi, 1 = 1,2, (49) 

for pseudoscalar, vector and scalar particles, respectively. Then, for pseudoscalar par- 
ticles (48) reads 



ML 2 (p) = E(p)Li(p)-j-^V(p-q L )(c-c--xs-s-)Li(q L ) , 



MLi(p) = ^(p)L 2 (p)-/ 7 ^\/(p-q)(c+c+-X<<)^(q) 



(50) 
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Here we have introduced as short-hand notation = s ±x , = c ±l . The quantity 
X = x(p, q) is given by (24). For vector particles system (48) takes the form 

o o r da 

M N« (p) = E(p) N? (p) + J j^V(p - q) 

•{(c pCq cr 6 - SpSq (6 ab X - 77V)) < (q) + V a c p c+ Ei (q)} , 

(51) 

MiV7(p) = E(p)^(p)+/^\/(p-q) 

■{(c+<r 6 - «(rV " ijV)) ^ (q) + ^c+c q & (q)} , 
where we have used the definitions 

V a = 77°(p, q) = &e?(p) , 77 a ee rf(p, q) = ^(q) , 5 ab ee ^(p, q) = e?(q)e?(p) • 
And for scalar particles one has 

ME 2 (p) = ^(p)Ei(p) 

+ / (0^ y (p - q){(* c P< - 4<) Si (q) +!z 6c p< (q)} > 

(52) 

MEi(p) = £(p)E 2 (p) 

+ / j0ys V (P ~ ^){(xc-c- - 8-s-) E 2 (q) +rfc;c- N b 2 (q)} . 

The relations (50)-(52) for bound state wave functions have a very compact form. In 
principle, they may be solved in dependence on the concrete form of the underlying 
potential. Thereby, in general an exact solution would be possible only by numerical 
calculations. But one can investigate also approximate solutions for different limiting 
procedures as, for example, nonrelativistic limits and the heavy quark mass limit. Let 
us add, that in [7] similar equations have been obtained for the oscillator potential, 
and the light meson mass spectrum has been calculated numerically. 
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5. Meson decay constants 



In the remaining sections of this paper we want to apply the bilocal meson model 
(9) with the relativistic covariant written kernel (11) to the investigation of heavy 
meson properties. Therefore we have to include into our QCD-motivated model the 
weak interaction. We will restrict ourselves here to the discussion of semileptonic weak 
processes. This allows us first of all to determine meson decay constants. Let us 
consider the quadratic part 

W§ = -^Tr(G s M) 2 (53) 

of the effective action (9). First of all we expand the bilocal field M. over creation (a^) 
and annihilation (an) operators 

= ^I^^I^^i^"^^ (54) 

+e- vs ¥a„(< i \V)tuHV)} • 

Here the sum runs over the set of quantum numbers H of hadrons contributing in the 
bilocal fields A4(x, y) . The bound state has the total momentum V = {ouh, V} , the 

energy ujh{V) = \pP + Mjj and the mass M H . The bound state vertex functions T H 
and Yh satisfy the Bethe-Salpeter equation (10) with kernel (11). 

Furthermore, we have to include the weak interactions into the effective action (53). 
The effective Lagrangian of semileptonic weak interaction has the form 

G F 

£semi = -^{VijiQiO^qj)^ + /i.e.} (55) 

with the leptonic current 

In = TO pin , I = e,n,T , M = 7^(1 + 75) , 

the elements Vij of the Kobayashi-Maskawa matrix and the Fermi constant Gf = 
1CT 5 / nip. Q denotes the column of (u, c, t) quarks and q - the column of (d, s, b) quarks. 
The lagrangian (55) can be incorporated into the bilocal action (53) by substituting 

M(x,y) -> M(x,y) + L(x,y) , (56) 

i.e., by adding to the bilocal field A4(x,y) the local weak leptonic current 

4(x, y) = ^5\x - y)Vje^^ , (57) 

where / = 0^1^ and Vl being the momentum of the leptonic pair. Then, the terms of 
interest standing in (53) after the substitution (56) and corresponding to semileptonic 



16 



weak interaction are 

W£u = -iN c Tr(GMGL) 

= —iN c J dxdydzdttr Gi'(t — x)A4i'k(x,y)Gji(y — z)Lj/k(z,t) 

Here the trace runs over Dirac and flavour indices. 

Now we are able to derive a formula for pseudoscalar meson decay constants. The 
matrix element for a decay of a meson ~ (QiQj) into a leptonic pair reads 

< lu\W^ mt \H tJ > = - t N c (27r) 4 t5^(V H -V L )^<lu\l fl \0> 

" / (^ tr ^{°^( g -^) f(q| ^ ) ^( 9 + ^)} • (58) 

Using the relations (33) for Green's functions and the definition of a bound state wave 
function (from (32) and (38)) 

. d( lo^ 7\fv_i<r>\^. /_ , V, 



we can write 



:/ ^G i (q-^)r(c i \V)G j ( q + ^) = ii> Hij (c i \V) , 



I ^^°^ Gi ( 9 ~^) mvH)Gi ( q+VB ^ 



(2tt) 4 7 l ^ l V 2 J v " n ' J V 2 J 
= ~ % \ J^^P^hM'Ph)} > (59) 

whereby ^(qj"^) = \l/(q| — Vh) ■ 111 analogy to (40) we introduce now the "dressed" 
o 

wave function \& and expand it in correspondence with (47) over the Lorentz matrices. 
Because we are interested in pseudoscalar mesons only, we take J = 1 and have then 
in the moving frame (using notation (49)) 

I (q| - V H ) = (ii (q| - V H ) - ff L 2 (q| - V H )^ . (60) 

Remember, that rf = V H /^Jv^ with V H = M H . Inserting (60) into (59) one gets after 
the calculation of the trace 

i^Vn J (cosi/i(q) • cosi/j(q) - sim/^q) ■ sinz/^q)) (L 2 )fly (q\V H ) , 

where z/(q) is defined in (20). Then, this expression is inserted into (57). Comparing 
the result with the general formula 

< lu\W^ mi \H l3 >= (2n) 4 i5 4 (V H -Vl^Fh^Vh < *|0 > 

one obtains for the decay constant of a pseudoscalar meson Hij at rest 

F Hij = ^ / 7^(^)^(q)cos(^(q) + ^(q)) . (61) 
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Notice, that in the rest frame the functions Li and L2 satisfy the system of equations 
(50). So we have for B and D mesons the relations 



AN C r dq ,0 



F D = T7~ I 77TZ^(^2) D cos(z/ c + v u 



4iV c f dq .0 
AN C f dq ,0 
AN C f dq ,0 



6. The relation between bilocal field approach and heavy 
quark effective theory 

For the description of the properties of heavy quarkonia one can employ the heavy 
quark effective theory, which has been developed recently [8 — 10]. Let us here consider 
mesons (Qq) consisting of a heavy quark Q and a light antiquark q. Now we make use 
of the fact that the quarkonium velocity is determined more or less by the velocity 
Vq of its heavy quark constituent Q. Then, the momentum of a bound state of velocity 
v M is given by 

V» = M/ , 

where M is the bound state mass. The latter is approximately equal to the heavy 
quark mass ijiq, i.e. 

M m rriQ . 

Then one can write for the heavy quark momentum 

Pq = m Q v Q = T*^ Pq = m Q v» + k» (62) 

with 

^ = (M - m Q X - p£ , (63) 

where p q is the light antiquark momentum. From (61) it follows that the heavy quark 
velocity is given by 

v% = v" + — k" . (64) 
Q m Q 

Therefore, in the limit itlq — > 00 one has Vq — > v^. In (62) and (63) the light antiquark 
momentum p^ is small as compared with V^ 1 , and it has only little influence on the 
direction of the heavy quark motion. As result one obtains that a meson bound state 
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containing a heavy quark can be considered by a theory describing heavy quark motion, 
for which the QCD corrections are small. Therefore the bound state is formed by a 
colour Coulomb field. 

Let us now assume that the bilocal field M. (x, y) contains a heavy quark (b or c) 
and a light antiquark (u, dors). Now, the integral kernel K, v (x — y), eq. (11), is defined 
with the help of the vector rj^, which has been introduced for two reasons. Firstly, it 
defines the transversality of the interaction given effectively by the phenomenological 
potential V^z^). And secondly, the vector rj^ determines the motion of the interaction 
potential together with the motion of the (Qq) bound state. In the limit of the heavy 
quark effective theory one has 

_ V, 

where v 2 — 1. Therefore the interaction kernel (11) takes in this limit the form 

K, v (x -y) = fViz^^iv ' Z )®i> , 

where z]| —v^v-z). We conclude that for considering the heavy quark effective theory 
one has to investigate the case of moving bound states in our bilocal field theory for 
nriQ — > oo. Concering the problems that have been discussed in this paper so far it was 
sufficient to consider the bound states at rest. Let us note the main modifications that 
appear if one works in an arbitrary reference frame. Instead of the 3-momentum q one 
has now q 1 , and q° is substituted by g". The product qcii changes to — ^ ± = — <j^7 M . 

Instead of |q| one has q 1 - = yj (g- 1 ) 2 , so that q 1 - = g^7^, q^ = q^/q ± - Then, for instance, 
eqs. (16) and (22) for S read 

E(^) = A{q L )q L + B(q^ = f + E(q ± )S~ 2 (q ± ) . 

For illustration let us now consider the Schwinger-Dyson equation in the heavy 
mass limit. Thereby, we will restrict ourselves to the investigation of the case (27), in 
which S(pg) = ttiq. Then, for S ±2 (q ± ) one obtains according to (20) 

S (q ) = — ±q 



E(q^) * E(q^) " 
The Schwinger-Dyson equation (15) for a heavy quark reads 

f d 4 q 

m Q = m° - itr J j^yV^fG^pQ - q)j , (65) 

where the trace runs over colour and Dirac indices. In the heavy mass limit we take 
inside the loop integral of (65) for the quark propagator the expression 



m Q {l + j)) - 4 
2rriQv(k — q) + q 2 + ie 
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Here we have used equation (62) and from (64) the relation k^/rriQ <C Vq ~ . 
Inserting (66) into (65) one gets 

m Q = m°-UN e l 'Av^) 7r ^ 5 

y i (2tt) 4 W J 2m Q v{k - q) + q 2 + ie 



m° + 2iV c / fV(^) , mQ 

J (2tt) 3 ^m^ 2 - 2m Q vk + (?) 



7. Semileptonic decays of heavy quarkonia in the limit of 
heavy quark effective theory 

Let us now consider semileptonic decays of heavy mesons. Therefore we have to 
investigate the cubic part of the effective action (9), thereby substituting one of the 
bilocal fields M by the local leptonic current L from (57): 

W£l = ^Tr[(G s M) 2 (G s L)] 

L(x4, x 5 )Gx(x 5 - x 6 )M(x 6 , xi)] . (67) 

Here the arguments in the integrand have been introduced according to fig.l. After 

rewriting (67) in momentum space and using the decomposition (54) for AA(x,y) one 

obtains for the matrix element describing a semileptonic decay of a B-meson H ib with 

a meson if'- in the final state 
j £ 

<(lv l )H' ii \W^ i \H ib > = -i^(2it)%V-V'-V L ) 



3 7 2(27r) 3 vW 
G_ 

V2' 



< lv\l tx \Q > .—V jb I» ji {V,V') (68) 



with 



.r^(k- yimW^Ck-liP) , (69) 



M = 7^(1 +75) and the latain indices 6, j — c,u and i — u,d,s,c indicating the quark 
content of the quantities. The momenta in the integral I b ^ have been introduced 
according to fig. 2. In (69) G b , Gj and Gi are the Green's functions of the quarks. For 
example, for j = c and % = u,d,s,c the matrix element (68) describes the decays 
B~ -> D°(lui), B° d -> D+(Zi4), 5° -> £>+(K>, and B~ -> (cc)(/z/,), respectively. 

For definiteness let us investigate the matrix element (68) for a semileptonic B- 
decay into a D-meson (j = c) . The calculation will be done as follows. Knowing from 



20 



sect. 6 the relation between bilocal field approach and heavy quark effective theory we 
will work at first in the rest frame. Only at the end we shall turn to the moving frame 
by substituting 70 by j] what corresponds for j] = f to the heavy quark mass limit. 

Proceeding in this way let us start with rewriting the integral (69). At first we 
insert the expressions (33) for the Green's functions into the latter: 

J i 2lT ) \ko- k\' + ie k - k 2 ' -ie) 

t Aj c) (k) | A<_ c) (k) \ r 
\ko — k[ c ^ + ie ko — k^ — ie) 

( Ag(k) Ag)(k) y Ahh {7Q) 

V k — k^ + ieko — k^ — ie) 

with 

kf/ 2 = M±E b (k) , 
k{% = M'± E c (k) , 
k[J 2 = ±£'j(k) . 

Now we represent the integrand of I^i as a sum of eight terms the numerators of which 
have the form 

tr[Ai fe) (k)7oO^Ai c) (k)7or„(k)7oA^(k)f^(k)] 

= tr^k^S-^kXi A±) Th (k)(± A±) K (k)(± A±)] . 



Here the last line is obtained by using (18) and (43), whereby r fulfils the Bethe- 
Salpeter equation (44). Because of the relations 

000 000 °u- 

A±ThA±= , A±ThA t =TI± t 

(cf. (42)) only two numerators remain. Then (70) is given by 

d 4 k 



o 

n*; (k) r H (k) a- 



.(k - kf - ie){ko - k[ c) + ie){ko - kf - ie) 
Uf_ (k) T H (k) A + 



{ko - k? + ie) (k - ki c) - ie) (k - k? + ie) , 
Now we can perform the k integration. The result reads 

dk 



IL(M,M') = itrf j^y 3 S^(k)0»S;\k 
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o 

u H ' (k) r H (k) a- 



[E c + E b -(M- M')](E C + Ei + M>) 



n?_ (k) h (k) a + ^ 



[E c + E b + (M- M')}(E C + Ei-M')J ' 

Here the expression in the brackets can be rewritten according to (45) as follows 

oo o ° , o oo , s o 2 o 

A-gg (k) A + r H (k) A_ _ A + fr g , (k) A_r g (k) A+ 

E c + + M ' £ c + Ei - M' 

L$h, (k) n"_ (k) £ + jr H , (k) n g + (k) 

E c + Ei + M> E c + Ei-M' 

Using (45) again and inserting the result in (71) one obtains 

dk 



1 bci 



(M, M') = itr | —^HkjOA-'l^tAiH' (k) A + * ff (k) A. 



o o 



- A + ^h> (k) A_#h (k) A+) . (72) 

In the following we will further rewrite the integral (72) in the heavy quark limit. 
To do this we have first of all to turn to the moving frame. We introduce the momenta 
as shown in fig. 2. In the moving frame we consider the decomposition fc M = + k^ 
of the momentum k^ = (ko, k) with respect to the momentum V of the initial boson. 
Furthermore, 70 has to be replaced by j] = f for the incoming b quark and by ify' = f 
for the outgoing c quark. Therefore in the moving frame the integral (72) takes the 
form 

rt 3 h 

lLAv,v 



■(l + f)^ H/ (k ± )i H (k ± )(l + ^) . (73) 

In this relation we have also made use of the fact that for heavy quarks the change to 
the moving frame leads immediately to the heavy quark effective theory. So for heavy 
quark fields Q with momentum V 1 = mv^ + k^ <C M one has 



Therefore the projector on the antipartilce vanishes: A^— ► 0. This means that for 



heavy quarks the only relevant contribution comes from A+^— > A+^ = (1 + 

Furthermore, we have used in (73) for the light quark q in the heavy quark limit 




A^-> 1/2. 
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For definiteness let us now consider the semileptonic decay of a pseudoscalar B 
meson into a heavy pseudoscalar meson of the type (ci), i — u,d, s, c. In this case we 
can use the decomposition (47) of the meson wave function which for moving fields 
takes the form 

h(k^\V) = 75 (£? + ^f)(A:». 
Inserting these decompositions into (73) one gets 



d 3 k A 

(2^ 

-(l + f + j + fflWft^v, v') , (74) 



W(h^\v,v') = L? (k^\V)L*{k^\V)+L? (k^V'^ik^V) 

+L»'(k ± \V , )L?(k ± \V) + L? (k L \V')L% (k^V) . (75) 

According to (20) we have 

Sr\k L )=dt-k L ^ i 

with 

c\ = cosiy i (k ± ) , = sinu i (k ± ) . 
After calculating the trace in (74) one obtains 

€ S) >,V) = i^(v + v'rj j^WdZ-XftW^M . (76) 

Inserting expression (76) for J^ 5)At into (68) one obtains for the matrix element in the 
case of a semileptonic decay of a pseudoscalar B meson into a pseudoscalar meson of 
the type (ci), i — u,d,s,c within the heavy quark effective theory the result 

<(lv)H' 3l \W^ mi \H lb > = ^(27ry5(V-V'-V L )<lu\l,\0> 

(v ■ v')(v + v'Y + £-(v • v')(v — v'Y] (77) 
with the Isgur-Wise functions [8] of the form 



t + (V ■ V = . ; 7=Vib 

v / MM 7 (27r) 2 v / ^ 7 v / 2 1 



T^ic^-s^Wik^v') , (78) 



(27T) 



Z-(v-v') = 0, 
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where W(k ± \v, v') is defined by (75). 
8. Summary and conclusions 

In this paper we have presented the bilocal field approach for relativistic covariant 
potential models in QCD. The main issue consists in the derivation of general integral 
expressions for meson properties as decay constants and semileptonic decay amplitudes. 
Therefore, it has been necessary to handle moving bound states. Doing this a special 
feature of our model has been important, namely, that the potential kernel (11) moves 
together with the bound state because of the presence of the vector 77^. Furthermore, 
to obtain semileptonic decay amplitudes we have established the relation between the 
bilocal field approach for our model and heavy quark effective theory. In this way we 
have been able to obtain an integral expression for the Isgur-Wise function. 

The formulas for the meson decay constants (61) and the Isgur-Wise function (78) 
appearing in the semileptonic decay amplitude (77) depend on the concrete form of the 
potential and contain trigonometric functions and meson wave functions. Thereby, the 
trigonometric functions fulfil together with the energy function the system of equations 
(23) which has been derived from the Schwinger-Dyson equation. For constituent 
quark masses one has the relations (27) and the system (23) simplifies significantly. 
The meson wave functions satisfy systems of equations (50)- (52) resulting from the 
Bethe-Salpeter equation. To obtain numerical results for the physical quantities one 
has to solve these systems of equations for concrete potentials. Work in this direction 
by using different approximations is in progress. 

We hope that the hadronization scheme with the use of the Foldy-Wouthuysen 
transform can be applied to the hadronization of quarkonia of the type (QQ) as well. 
This would allow one to describe also nonleptonic decays of heavy mesons in the same 
fashion. 
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Tables 



Potential 


V(r) 


V( P ) 


Coulomb 


4 a s 
3 r 




Linear 


ar 




Oscillator 


br 2 


-6(2tt) 3 A p 5( P ) 


Yukawa 


r 




NJL 


V 5(r) 




Constant 


V 


K,(2tt) 3 5( P ) 



Tab.l: Interaction potentials in x space and momentum space. 
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Figure captions 



Fig.l: Diagram for the semileptonic decay of a heavy meson A4(x2,xs) into a heavy 
meson AA{x§,Xi) and a leptonic current L(x±, x 5 ) in bilocal field theory. 



Fig. 2: Diagram corresponding to the integral J^, eq.(69), figurating in the semileptonic 
decay amplitude (68). 
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